Discovery of a possibly single blue supergiant star in the 
intra-cluster region of Virgo cluster of galaxies 
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ABSTRACT 

IC 3418 is a dwarf irregular galaxy falling into the Virgo cluster, and a 17 kpc long trail is 
seen behind the galaxy, which is considered to have formed due to ram pressure stripping. The 
trail contains compact knots and diffuse blobs of ultraviolet and blue optical emission and, thus, 
it is a clear site of recent star formation but in an unusual environment, surrounded by a million 
degree intra-cluster medium. We report on our optical spectroscopy of a compact source in the 
trail, SDSS J122952. 66+112227. 8, and show that the optical spectrum is dominated by emission 
from a massive blue supergiant star. If confirmed, our report would mark the farthest star with 
spectroscopic observation. We interpret that a massive 0-type star formed in situ in the trail 
has evolved recently out of the main sequence into this blue supergiant phase, and now lacks any 
detectable spectral sign of its associated H ii region. We argue that turbulence within the ram 
pressure striped gaseous trail may play a dominant role for the star formation within such trails. 



Subject headings: galaxies: clusters: individual (Virgo) 
medium — stars: formation — supergiants 

Introduction 



As galaxies fall into the huge gravitational po- 
tential of clusters of galaxies with high velocity, 
interaction with the hot intra-cluster medium 
strips off cool gas from the main body of the 
i nfalling galaxies into the intra-cluste r region 
(JGunn fc Gott Ill972l : lwlmer et al.ll200l[) . Study- 
ing fate of this stripped gas, ionization in contact 
with million degree hot intra-cluster medium, or 
condensation to form new stars, is of great in- 
terest. The Virgo cluster of galaxies (hereafter, 
the Virgo cluster) is t he nearest mas sive cluster 
of galaxies (16.5 Mpc; iMei et al.l 120071 ) . and is an 
ideal laboratory to investigate details of such pro- 
cesses. We focus on IC 3418, a galaxy likely falling 
into the Virgo cluster for the first time at a very 
high speed (nearly 1000 km s^^; IVollmer et al 



■ galaxies: individual (IC 3418) — intergalactic 



20051 ) UV continu um and Ha emission line im- 
ages ([Hester et al 



20101: iFumagain et"aLl 1201 1[) 



Both of these emissions are clear signs of recent 
star formation, and therefore, the trail is the best 
target to explore details of young massive stel- 
lar population in the intra-cluster region. Stellar 
spectroscopy is unique and useful not only because 
it provides us basic properties of the stellar prop- 
erties (e.g., age, metallicity, luminosity) but also 
because we can explore a potentially new mode 
of star formation in this special condition. We 
expect that such star formation is quite different 
from those in Milky Way, and such study will even- 
tually contribute to understanding of galaxy evo- 



lution under cluster environment (Gunn fc Gott 



20011 ) ■ especially on its 17 kpc long trail which 
is considered to have formed behind the galax^ 



due to ram pressure s tripping (jHester et al 



aiaxy 
2Qlo[ 



19721 ) and intra-cluster stellar population. 

In this letter, we report on our spectroscopic 
studies of SDSS J122952.66-M12227.8 (hereafter, 
SDSSJ1229-I-1122), a compact source of optical 
emission within the star-forming trail of IC 3418. 



Fumagalli et al.ll2011 ). The trail is m anifested as a 
chain of blobs seen both in GALEX (Martin et al. 



2. Spectroscopic Observation and Pho- 
tometry on Archival Images 

FOC AS spectrograp h at Subaru telescope 
(JKashikawa et al.l 120021 ) was used, on the night 
of 27 April 2011, in multi-slit spectroscopy mode 
to study the nature of the trail of IC 3418. One 
of our slitlets was put on a faint {g'j^^—23.0 
mag) blue compact object catalogued earlier as 
SDSSJ1229+1122 and located near the far end 
of the trail (Figures 1 and 2). Low-resolution 
spectrum was taken with 0.8"-wide slit and 
"300B" grism for spectral resolution R{= X/SX) 
of ~ 700 to cover 4000A-7000A, and the total 
exposure time was 6000 s. Galactic extinction for 
E{B — y)=0.03 magnitude was corrected based 
on the NASA Extragalactic Database (NED) in- 
formation. Absolute flux scale of the spectrum 
was re-calibrated based on our photometric data 
(see below), and a small slit- loss correction (by 
7%) was applied. Additionally, a shorter exposure 
(3600 s) but with relatively higher ("medium") 
spectral resolution {R ~ 3750) spectrum was ob- 
tained with the same slit width and "VPH680" 
grism, resulting in full width at half maximum 
(hereafter, FWHM) velocity resolution of ~ 80 
km s^^. Low-resolution spectrum is shown in Fig- 
ure 3, and measured spectroscopic properties on 
both spectra are summarized in Table 1. 

We also performed photometric analysis on 
archival images to supplement our spectroscopy 
(Figures 1 and 2). Canada- France-Hawaii Tele- 



scope (CFHT) Megacam dBoulade et al.l 



2003 ) 



20081 ) 



images were taken from Megapipe ( Gywnl 
at g' , i' , and z' bands, and another image at u* 
band was locally stacked based on individually 
pre-calibrated archival images available at CFHT 
Science Data Archive. All those images were orig- 
inally taken for New Ge neration Virgo cluste r 
Survey (NGVS) project (JFerrarese et al.l 12012"). 
SDSSJ1229+1122 has a faint nearby source within 
1" from our spectroscopic target source, which was 
not resolved in SDSS image. SDSSJ1229-^1122 
was unresolved in both the i' and z' band images 
having the best FWHM seeing size of 0.7" (~ 60 pc 
at the distance of the Virgo cluster) . Multi-color 
photometry was performed through double Gaus- 
sian fitting on the images to properly remove any 
contamination from the nearby source, and the 
same amount of Galactic extinction is corrected 



as for the spectra. In addition. Ha image was ob- 
tained with EF0SC2 attached to NTT telescope 
on 2011 April 29. Narrow "Ha#692" fiher cov- 
ering Ha within velocity range of —780 to 2060 
km s~^ and broad-band "R#642" filter are used 
for stellar continuum subtracted Ha line image. 
Total exposure times were 1800 s and 600 s for Ha 
and R band images, respectively. Flux calibration 
was performed with a spectroscopic standard star 
LTT3218. Our point source detection limit (3(j) is 
1.8x10"^^ erg s"-"^ cm~^. Note that this is an inde- 



pendent im age from that pre s ented in lHester et al.l 
(|2010l ) and iFumagalli et al.l (|201l[ ). Also, deeper 
GALEX NUV image (1750A-2800A) was created 
by combining the GA LEX Ultraviolet Vir go Clus- 
ter Survey (GUViCS) (iBoselli et al.ll201ll) a nd the 
Deep Imaging Survey CDIS) (JMartin et al.llioosi) . 
The total integration time is 20672 s. All the im- 
ages are shown in Figures 1 and 2, and optical 
photometry data (Table 1) are overlaid on spectra 
(Figures 3 and 4). 

3. Results 

The low-resolution spectrum of SDSSJ1229-I-1122 
can be characterized by strong Ha (at equiva- 
lent width (EW) of — 34A) and weak H/3 emission 
lines on blue continuum. The Ha line is slightly 
resolved with our medium-resolution spectrum, 
having FWHM of '^ 164 km^^ (corrected for the 
instrumental resolution), but further line profile 
analysis was not possible. Our spectroscopic Ha 
flux (5.0 X 10^^^ erg s^^ cm^'^) is below detec- 
tion limit of our Ha image (Figure 2). The Ha 
flux is < 10 times fainter th an nearby compact 
Ha emitters within t he trail ( Hester et al.ll2010t 
Fumagalli et al.l 1201 It see also Figure 2). Our 
medium-resolution spectrum shows that the he- 



liocentric velocity is —99 km s 
close to that of IC 3418 (38 km s~ 



which is very 



Gavazzi et al 



Their similar recession velocities as well as 
spatial association of SDSSJ1229-H1122 with the 
trail of IC 3418 (Figure 1) strongly support that 
they are physically associated to each other within 
the Virgo cluster. At a distance o f 16.5 Mpc (or 



m 



M = 31.1; iMei et al.l [20071 ). the absolute 



magnitude (Afy) is ~ —8.3 and Ha luminosity is 
1.6 X 10"^^ erg s~^. The Ha luminosity is about 5 
- 100 times fainter than known intra-cluster H 11 



regions around NG C 4388 () Gerhard et al.M2002 : 



Cortese et al.l 120041 ) and about 10 times fainter 



than ones around VCC 1249 / M 49 merging system 
(JArrigo ni Battaia et al.ll2012l ) in the Virgo cluster. 
In one of these known H ii regions, where detailed 
spectroscopic analyses were made, only a few O- 
type stars were estimated t o be powerful e nough 
for the observed luminosity ( Gerhard et al.l l2002'). 

We find that the emission-line ratios of SDSSJ1229- 
are unusual for photo- or shock-ionized nebulae. 
Other than Ha and H/3, no nebular forbidden lines 
were detected, including ones that are typically 
bright in those nebulae ([O in] A 5007A, [N ii]AA 
6548,65831, and [S ii]AA 6717,6731A) (Figure 3). 
Their upper limits were measured assuming the 
same line widths and redshift of Balmer emis- 
sion lines (Table 1). Stringent upper limits on 
low ionization [N ii] and [S ii] and undetected 
higher ionization [O in] lines strongly indicate that 
the emission originates nei ther from H n regions 
nor s hock-ionized nebulae (|Veilleux fc Osterbrock 
19871 ) but from hot stellar atmosphere. We also 
note that other emission lines seen typically in Lu- 
minous Blue Variable (LBV) stars besides Balmer 
lines, such as [He i]A5878A and [Fe II]A5018A, 
were not detected. 

Such a luminous compact object with blue con- 
tinuum should either be a young stellar associa- 
tion or a single blue (super)giant s tar. We first 
searched a stellar spectrum library (|Pickleslll998l ) 
for matching the continuum (below Ha and ex- 
cluding H/3 emission lines) and the broad-band 
photometric data, and found that A2 supergiant 
(I) is the best match to the observation (Fig- 
ure 4). Similar blue continuum can be repro- 
duced by A5 dwarf (V), but its fainter expected 
u* fiux than the observation as well as intrin- 
sically faint luminosity does not agree with the 
observation (Figure 4). Then we searched star- 
burst99 models of in stantaneous star formation 
( Leitherer et al.lll999l l for the same purpose. Con- 
tinuous star formation models are not considered 
since they should always show H n-region like 
spectrum, unlike our observation. We chose a 
model of 0.4 solar metallicity (Z = 0.008) based on 
metallically measurement with oth er intra-cluster 



H II regions in the Virgo cluster (jGerhard et al. 



20021 : lArrigoni Battaia et al.ll2012[ ). Shape of ini- 
tial mass function (IMF) can not be examined, 
and we only describe here cases with IMF slope 
(a)=2.35 and Afuppcr = 100 Mq since basic re- 
sults hold with other IMF parameter sets. Al- 



though the model grid provided bv lLeitherer et al 
(|l999l ) is not fine enough to find statistically satis- 
factory models for all photometry data, we found 
three possible burst ages to match the observa- 
tion: 6 Myr, 16 Myr, and 20 Myr (Figure 4). 
We found three solutions because, in addition to 
-1122iiain sequence stars that become monotonically 
redder with age, enhanced red supergiant contri- 
bution quickly makes the color redder just around 
~ 10 Myr (or between 6-16 Myr in the assumed 
model parameters) when most 0-type stars evolve 
away from the main sequence. Note that the very 
similar trend of color change is seen for all different 
IMF parameter sets. In summary, both possibili- 
ties of a blue supergiant (hereafter, BSG) star and 
young (6-20 Myr) stellar association are equally 
good to reproduce the observed blue continuum 
spectrum and broad-band photometry. 

4. Discussions 



4.1. SDSSJ1229+1122: A compact young 
stellar association or a single blue su- 
pergiant star? 

In the model of instantaneous starburst, the 
emission line component of SDSSJ1229-1-1122 
should originate from nebula ionized by young 
massive stars within the stellar association. Since 
EW(Ha) monotonically decreases with the burst 
age, only the ~ 6 Myr model among the three 
possibilities can reproduce rather large observed 
EW(Ha) (— 3 4A) due to rema ining main sequence 
0-type stars ( Leitherer et al.lll999 '). The contin- 
uum luminosity then requires stellar mass of '^ 10** 
Mq, and expected number of 0-type stars there 
is several - several 10s depending on choices of 
the IMF parameters, producing luminous H ii re- 
gion spectrum that is not observed. Next pos- 
sibility is that a stellar association of ~ 16- 
20 Myr dominates the observed continuum and 
extra Ha emitting source(s) in addition to regu- 
lar population produce strong Ha. Possible extra 
emission-line components are classic al Be stars 
( Lamers et al.lll998l : IPaul et al.ll201 Z). supergiant 
Be stars, and supergia,nt A stars (iTully &: Wolfil 
1984 iMcCarthvet al.lll997l ). Ah of them show 



prominent Ha emission without nebular lines, and 
their EW(Ha) alo ne is already comparable to the 



observed one fe.g.. lTullv fc WoM1984HPaul et al 



20121 ) , indicating that those stars dominate the ob- 



served spectrum since underlying spectrum from 
regular stellar association, if any, will dilute the 
observed width. The observed absolute magnitude 
{My ^ —8.3) requires either a singl e BSG or 50- 
200 c l assical Be stars ( class III-V) ( Lamers et al.l 
llQQSt IPaul et al.l 120121 ) since they are fainter by 
4-6 magnitude than BSGs. If SDSSJ1229+1122 
is made entirely of classical Be stars, such larger 
number of short-lived objects within a compact 
(^ 60 PC across) stellar association requires a 
very unlikely star formation episode. Instead, if 
it is a BSG, both the absolute magnitude and 
EW(I Ia) are within the range of AO I stars in 
SMC (JTullv &: WolfJll984l) although the observed 
EW(H/3) is larger by a factor of ~ 2. Since the 
EW in BSGs is a funct ion of metallicity, l umi- 
nosity, and stehar type (JTullv fc WolfJll984l) . we 
believe that our observation is within the probable 
range of BSGs given uncertainties of detailed stel- 
lar parameters for SDSSJ1229-hll22. Observed 
intrinsically broad (FWHM of ~ 164 km"^) Ha 
emission hints presence of strong stellar wind 
and/or fast rotating equatorial disk, although 
their char acteristic profile (e.g., P-Cygni type 



with UV-bright, Ha emitting, an d/or optically 



for BSGs; iMcCarthv et al.l 119971: iPrvbilla et al. 



20061 ) could not be discerned. Therefore a model 



of single BSG as the optically dominant sources 
within SDSSJ1229-I-1122 is strongly preferred over 
a model of young stellar association. Note, how- 
ever, that we can not distinguish between case of a 
single star or collection of a few such stars based on 
our observation, although only a single such star is 
bright enough. We expect that SDSSJ1229-I-1122 
was an H ii region about a few - several Myr 
ago and similar in nature to other known H ii 
regions in the trail of IC 3418 ( Hester et al.luOlOt 
Fumagalli et al.l 2011 ) and ones within the Virgo 



cluster (iGerhard et al.ll2002t ICortese et al.l 12004: 
Arrigoni Battaia et al.ll2012l) 



4.2. Star formation characteristics within 
the cluster environment 

Stars usually form in dusty giant molecu- 
lar gas clouds, but such clouds are too tightly 
bound with the galaxy to be stripped off by ram 
pressure stripping and thrown into the gaseous 



trail of cluster-infalling galaxies ([Hester et al 
2010l ). However, a few rare trails including 



blue young star-forming clumps (ICortese et al 
20071: IV oshida et al. 2008; Hes ter et all l20ld 



Fumagalh et al.|201lHArrigoni Battaia et al.ll2012l) 



Fumagalli et al.l (|201ll ) analyzed stellar proper- 
ties of IC 3418 based on optical and GALEX- 
UV spectral energy distribution (SED), and ar- 
gue that the star formation there has been trun- 
cated about 200 Myr ago due to ram pressure 
stripping. As massive stars are short-lived, the 
0-type star, which has evolved into a BSG by 
now, must have formed a few 10 Myr ago within 
SDSSJ1229-t-1122, i.e., it formed after the trun- 
cation of star formation within IC 3418. This 
supports an idea that the star have formed in situ 
within the trail. 

SDSSJ1229-fll22 is located wit hin a diffuse 
GALE X NUV blob, labeled "D3" bv lHester et"al 
( 20101 ) or a "lower surface-brightness and vi- 
sually more dif f use fi lament", labeled "Fl" by 



Fumagalh er ahl (l201ll). The blo b extends ~ 10 



or ~ 800 pc. iFumagalli et al.l (|2011l) reported 



one of IC 3418 and others that extend up to 
90 kpc from the parent galaxy are seen decorated 



2.3 magnitude brighter g' band magnitude for Fl 
than our measurement on SDSSJ1229-|~1122, and 
analyzed that the overall SED including GALEX- 
UV bands can be reproduced by young stellar 
population with age of 6401^^44 Myr and mass 
of 3.48 X 10^ Mq. However, the CFHT image 
shows only a few more fainter compact sources 
besides SDSSJ1229-hll22 there (Figure 2), indi- 
cating that there is an unresolved diffuse young 
stellar population within the blob. We speculate 
that there are much more unexplored population 
of blue stellar population within the intra-cluster 
region, not only along the trail of IC 3418 but also 
in other similar trails in the Virgo cluster. 

It is likely that star formation in such a spe- 
cial environment is quite different from typical 
cases in, e.g., our own Milky Way galaxy. Differ- 
ence in temperature and relative velocity between 
the star-forming molecular clouds and the ambi- 
ent medium may be million degree and ~ 1000 
km s~^, respectively. Being far from the par- 
ent galaxy, especially at the far end of the trail, 
role of turbulence may be dominant over grav- 
ity. Turbulence in the trail may create eddies to 
form dense cloud clumps which would cool very 
fast and subsequently collapse under their own 
gravity to form stars. Although the model of 
turbulence-driven star formation seems promising 



(JHester et all l2010t iFumagalli efsd] l201ll) . con- 
straining detailed physical mechanism of the star 
formation and characterizing the turbulence re- 
quires further investigations of stellar population 
and cool molecular gas there. 

5. Implication for Future Extragalactic 
Stellar Astrophysical Studies 

Various kinds of intra-cluster stellar popula- 
tion have been known in the Virgo cluster: In 
additi on to diffuse intra- cluster optical continuum 
light (|Mihos et al.l 120051 ). i ndividual sources like 
intra-cluster H il regions ( Gerhard et al.l 2002t 
Cortese et al.ll2004l: lArrigoni Battaia et al.ll2012l). 
old ( > 10 Gyr) intra-cluster stars dWillianis et al 



2007 ) , and intra-cluster plan etary nebulae (jFeldmeier 



2004t lAguerri et al.l |2005[) are known. O 
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Ferrarese, L. et al., 2012, ApJS, 200, id. 4 

Fumagalli, M., Gavazzi, G., Scaramella, R., & 
ITCiadrSia^^^'tti, R 2011, A&A, 528, 46 



study, independent of the interpretation that 
SDSSJ1229-t-1122 contains a single or a few more 
stars, has added BSGs to be residing in the 
intra-cluster region. Previously, BSG candidates 
have been phot ometrically identified in M 100 
(JHill et al.lll998l ) in the Virgo cluster, and spec- 
troscopically identified in NGC 3621 at a distance 
of 6.7 Mpc (Bresohn et al. 2001). Being in the 
Virgo cluster, and if confirmed by future obser- 
vations to be a single star, SDSSJ1229+1122 will 
be the most distant star discovered from spec- 
troscopic observations. We demonstrated for the 
first time that stellar spectroscopy up to distance 
of the Virgo cluster is indeed f easible as antic- 
ipated nearly two decade ago (IKudritzki et al 



1995; 'McCar thy et al.lll997l: iBresolin et al.l 12001 



Kudritzki 20101 ). However, to reveal its true 



na- 
ture by detailed quantitative spectroscopic analy- 
ses, spectra with higher signal-to-noise ratio and 
higher resolution below H/3 (< 4800 A) as well as 
Ha are required. To exploit its full potential for 
accurate distance measurements and star forma- 
tion studies in the Virgo cluster and beyond, we do 
need the next-generation giant telescopes such as 
Thirty Meter Telescope (TMT) and/or European 
Extre mely Large Telescope (E-ELT) (jKudritzkJ 
l201fl[ ). 
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Table 1: Photometric and spectroscopic character- 
istics of SDSSJ1229+1122 



Megacam astrometry 



RA (2000) 
12'^29'"52^69 



Megacam photometry^ 

u* 

9' 



Bessel photometry^'^ 

V 



DEC (2000) 
+11°22^28.0" 

Corrected'^ 
23.15 ±0.02 
22.85 ±0.01 
23.03 ±0.02 
23.12 ±0.04 

Absolute^ 
Mv = -8.25 



Observed 
23.29 ±0.02 
22.96 ±0.01 
23.09 ±0.02 
23.16 ±0.04 

Apparent 
22.85 ±0.01 



FOCAS low-resolution spectroscopy 

Ha flux2 

Ha equivalent width 

H/3 flux2 

H/3 equivalent width 



(5.04 ± 0.08) X 10^1^ erg s^^ cm^^ 

-33.7 ±0.7 A 
(2.08 ± 0.22) xlO-" erg s"! cm^^ 

-6.7 ±0.7 A 



Emission line ratios 
[OIII]A5007/H/3 
[NH]A6583/Ha 
[SH]A6717+6731/Ha 



3cr upper limit 



<0.28 

< 0.028 

< 0.031 



FOCAS medium-resolution spectroscopy 




Ha heliocentric velocity 


-99 ±6 kms"^ 


Ha intrinsic line width 


164±22kms^i FWHM 



^In Megacam SDSS {AB) system. 

^Galactic extinction corrected for E(B — y)=0.03. 

^Bessel system in Vega unit, converted based on our spectrum. 

■^Distance modulus o{ m — M = 31.1 is used. 
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Fig. 1.— GALEX NUV (a) and CFHT g' band (b) images of IC 3418 and its trail including SDSSJ1229+1122 
(marked by two arrows) in linear intensity scale. An inset g' band image within (b) shows fainter tidal features 
of IC 3418 in square-root intensity scale. 
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Fig. 2. — Zoomed-in optical views around SDSSJ1229+1122. Upper panels show ^ V view around 
SDSSJ1229+1122 in u* , g' , i', and continuum-subtracted Ha line ((a) - (d)). Lower panels show the 
same images but just for central ^ 10" region ((e) - (h)). All images are shown in linear intensity scale. 
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Fig. 3. — Top panel: Low-resolution optical spectrum and u*, g', and i' band photometry [blue with error 
bars) of SDSSJ1229+1122. One-sigma noise spectrum is shown near the bottom of all plots, bottom left: 
The same spectrum of top panel but around H/3 and [O in] . Multi-Gaussian fitting (plus linear function for 
the continuum) results and the line wavelengths are overlaid with cyan broken lines. Neither [O iii]A4959 
nor A5007 are detected but their wavelengths expected for the redshift of H/3 are marked, bottom right: The 
same spectrum of top but around Ha, [N ii], and [S ll]. Multi- Gaussian fitting results are shown in the same 
way as for the bottom left. Neither [N ii]A6583, [S ii]A6717, nor A6731 are detected but their wavelengths 
expected for the redshift of Ha are marked. 
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Fig. 4. — Comparison of observed low-resolution spectrum and photometry data for u* , 5', i', and z' bands 
{blue with error bars) of SDSSJ1229+1122 [top) with model spectra. Model spectra are scaled at g' band to 
the observation, and then offset by —0.25 in log flux scale along with the photometry data for clarity of the 
figure. As labeled, second, third, and fourth spectra are starhurst99 instantaneous models of 6 Myr, 16 Myr, 
and 20 Myr, and following two spectra are of A2 I and A5 V stars, respectively. Expected u* fluxes for each 
model are marked with red crosses. 
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